Dynamic Coarse-Graining Approach to Quantum Field Theory 
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We build quantum field theory on the thermodynamic master equation for dissipative quantum 
systems. The vacuum is represented by a thermodynamic equilibrium state in the low-temperature 
limit. All regularization is consistently provided by a friction mechanism; with decreasing friction 
parameter, only degrees of freedom on shorter and shorter length scales are damped out of a quantum 
field theory. No divergent integrals need to be manipulated. Renormalization occurs as a tool to 
refine perturbation expansions, not to remove divergences. Relativistic covariance is recovered in 
the final results. We illustrate the proposed thermodynamic approach to quantum fields for the 
<p 4 theory by calculating the propagator and the f} function, and we offer some suggestions on its 
application to gauge theories. 



I. INTRODUCTION 

Renormalization plays a key role in quantum field the- 
ory. The renormalization group is intrinsically related 
to the successive elimination of degrees of freedom and 
should hence be expected to result in coarse-grained 
equations. Quantum field theory hence belongs into the 
world of multiscale modeling and coarse graining of dy- 
namic systems. The field idealization implies the need 
to bridge a wide range of length scales and, in relativis- 
tic quantum field theory, clearly also in time scales. We 
hence need a proper framework to implement dynamic 
renormalization-group theory. 

The proper setting for multiscale modeling and coarse 
graining is statistical nonequilibrium thermodynamics. 
Whereas, in the context of quantum field theory, one usu- 
ally does not think of irreversible processes, the creation 
and annihilation of particle-antiparticle pairs at inacces- 
sibly small length and time scales, for example, certainly 
is beyond our control and should hence be regarded and 
treated as an irreversible process. We hence propose to 
make use of the powerful machinery of nonequilibrium 
thermodynamics instead of inventing a set of sophisti- 
cated rules to eliminate various kinds of infinities from 
an approach designed for reversible systems. 

As a result of the thermodynamic approach, we avoid 
the possible concern that renormalization "is simply a 
way to sweep the difficulties of the divergences of [quan- 
tum] electrodynamics under the rug," as expressed by 
Feynman in a catchy metaphorical statement in his No- 
bel lecture (1965). Or, in the words of the insistent critic 
Dirac "the quantum mechanics that most physicists 
are using nowadays [in quantum field theory] is just a set 
of working rules, and not a complete dynamical theory 
at all." Whereas Dirac felt the need for a different type 
of Hamiltonian, we here suggest to address the intrinsic 
irreversibility associated with the field idealization and 
renormalization in an appropriate manner. Maybe such a 
thermodynamic approach introducing irreversibility into 
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quantum field theory could eventually provide "some re- 
ally drastic changes" in the equations, as demanded by 
Dirac (see pp. 36-37 of Q). 

We start by introducing some relevant background in- 
formation: The quantum master equation required to 
formulate the dissipative dynamics of quantum systems 
and the renormalization-group approach employed to re- 
fine the results of perturbation theory. We then develop 
the thermodynamic approach for a scalar quantum field 
with quartic interactions. It is shown how the interac- 
tions can be handled by a straightforward, thermody- 
namically consistent perturbation theory that takes de- 
tailed balance into account; first and second order re- 
sults for the propagator and a four-point correlation are 
calculated and discussed in great detail. The steps re- 
quired to generalize the dissipative approach from scalar 
to more complicated fields, such as Yang-Mills gauge 
fields, are sketched in a further section. We conclude with 
a brief summary and some further remarks on the coarse- 
graining approach to quantum field theory. Useful details 
are compiled in a number of appendices. Throughout this 
paper, we use units with K = c = 1 where H is the reduced 
Planck constant and c is the speed of light. 



II. QUANTUM MASTER EQUATION 

In this section, we compile the equations required 
to describe the mixed reversible-irreversible evolution of 
quantum systems so that we can later formulate a dissi- 
pative smoothing mechanism for quantum field theories. 
After introducing the general framework, we focus on 
the equations obtained by linearization around an equi- 
librium state, on detailed balance, and the resulting guid- 
ance for the development of perturbation theory. 



A. Thermodynamic master equation 

A standard approach to quantum dissipation is based 
on linear quantum master equations for density matri- 
ces 0, The form of the most general linear quan- 
tum master equation that is guaranteed to preserve the 
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trace and positive-semidefiniteness of any initial density 
matrix has been determined by Lindblad [5j. However, 
master equations of the Lindblad form are not applicable 
to arbitrarily low temperatures. As we wish to construct 
quantum field theories in the limit of zero temperature 
and vanishing dissipation rate, we need to pay attention 
to the thermodynamic consistency of the quantum mas- 
ter equation to render it applicable for low temperatures 
and weak friction jg|. 

Motivated by a failure of the "quantum regression 
hypothesis," Grabert revisited the standard projection- 
operator derivation of quantum master equations [7[ with 
a relevant density matrix of the exponential form, where 
the deviation from the Hamiltonian in the exponent can 
be interpreted as the thermodynamic force operator con- 
jugate to the density matrix. In the Markovian limit, 
the resulting equation for the evolution of the density 
matrix or statistical operator p is of the nonlinear form 
(see Eq. (5.22) of Q), 

3,k 

- k B T c ^2j jk [Qj, [Q k ,p]}- (1) 

The reversible first term on the right-hand side of Eq. (JTJ 
has the well-known representation in terms of the com- 
mutator and the Hamiltonian H. The additional irre- 
versible term is formulated in terms of a suitable matrix 
of damping coefficients jjk describing the strength of the 
dissipation and the observables Qj describing the inter- 
action between the quantum subsystem and its environ- 
ment which, in the case of quantum field theory, is a bath 
representing the unresolvable or eliminated small-scale 
features. As usual, k B is Boltzmann's constant and T e 
is the temperature. Nonlinearity arises from the super- 
operator 

K P A= [ p u A P 1 ~ u du, (2) 
Jo 

which adds a factor of p to the observable A with a proper 
treatment of ordering problems. Note that the temper- 
ature T e is the only parameter characterizing the state 
of the environment, which can hence be regarded as a 
heat bath. Equation (TT|) may be addressed as a ther- 
modynamic master equation because it has been derived 
with a relevant density matrix characterized in terms of 
a thermodynamic force operator and because, as a con- 
sequence, it is consistent with the fluctuation-dissipation 
theorem. 

The quantum master equation ([1]) and its generaliza- 
tion for the coupling to arbitrary classical nonequilibrium 
systems as environments has recently been obtained in 
9] . Starting point is the geometric formulation of classi- 
cal nonequilibrium thermodynamics in terms of Poisson 
and dissipative brackets [Iol - |T2| . The quantum general- 
ization is obtained by Dirac's method of classical analogy 



[13j |: Poisson brackets are replaced by commutators, dis- 
sipative brackets are replaced by canonical correlations 
of commutators. The appealing properties of the nonlin- 
ear quantum master equation (TTJ) have been discussed in 
[Til ]. As the purely phenomenological and formal projec- 
tion operator approaches to quantum dissipation lead to 
the same nonlinear thermodynamic master equation ap- 
plicable in the low-temperature regime, we have a save 
starting point for our dissipative approach to quantum 
field theory. 

B. Linearized master equation 

In order to clarify the nonlinear nature of the thermo- 
dynamic quantum master equation Jl} and to linearize it 
around its equilibrium solution 

Pcq = tr e~P H ' 

at inverse temperature j3 = (fceTe) -1 , we introduce the 
thermodynamic driving force operator, 

p = k B T c (Inp- lnpcq). (4) 

We can then write 

K p [Qk,H] =]C p [Qk,ij] - k B T e [fC p Qk, In p\- (5) 

With this simple consequence of the definition of the ther- 
modynamic driving force operator and the useful identity 

[A,p] = [K p A,\np]=1Cp\AMpl (6) 

which follows from looking at arbitrary matrix elements 
formed with the eigenstates of the density matrix and 
performing the elementary integration over u in Eq. , 
we can rewrite the quantum master equation (TTJ) in the 
more compact form 

f t =i[p,H]-Y,l3k[Q3^ P [Qk^]\- (7) 

In view of the logarithmic form of p given in Eq. ((4]), this 
equation is clearly nonlinear in p. It moreover justifies 
the interpretation of p as a thermodynamic driving force 
operator. 

In linearizing around equilibrium, the driving force p 
is small so that K, p in Eq. (J7]) can be evaluated with 
the equilibrium density matrix; in that case, we use the 
symbol K. From the linearized relationship between the 
deviation from the equilibrium density matrix and the 
thermodynamic driving force operator, 

p = exp{lnp cq + fifi} w p oq + K./3p, (8) 

we finally obtain the linearized thermodynamic master 
equation in the form 

^ = i[p,H] - k B T e Y,'YjklQ j ,)C[Qk,)C- l p}}- (9) 
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Note that the master equation (9) is not of the Lindblad 
form 5]. For large deviations from equilibrium, positive- 
semidefmiteness of the density matrix is only guaranteed 
by the nonlinear terms. For small deviations from equi- 
librium, we can safely use Eq. . 



C. Adjointness properties 

For the general development of the theory of linearized 
thermodynamic quantum master equations, it is conve- 
nient to look at the time evolution of the observable 
p = K,~ 1 p. In other words, p is the relative deviation of 
the density matrix p from its equilibrium form, p cq . By 
operating with /C _1 on the linearized quantum master 
equation © , the evolution equation for the observable p 
is given by 

^ = Cp = i[p,H] - k B T c ^2"( jk K.- 1 [Q j ,K.[Q k ,p}}. 

j,k 

(10) 

Of course, p = 1 is found to be the steady state solution 
of this equation. 

The average of any observable A can be obtained as 

tr(Ap)=tx(A)Cp) = (A;p), (11) 

where the second identity is the definition of the equi- 
librium canonical correlation (see Eq. (4.1.12) of [15j or 
Eq. (3) of @). The time evolution of an average can be 
obtained from the master equation {9]) for the evolution 
of p or from Eq. (JTUJ) for p. According to the property 

(A-£B)=(CA;B), (12) 

where the adjoint operator L follows from an elementary 
calculation, 

LA = -i[A, H] - k B T c ]T 7kj£~ 1 [Qj>fc[Qk,A]], (13) 

there exists a third possibility to obtain the evolution 
of an average: one can use the time-dependent observ- 
able A obtained from dA/dt = LA and a constant initial 
density matrix p or p in Eq. (|11[) . The reversible con- 
tribution to the evolution equation for A corresponds to 
the usual Heisenberg equation for observables. The op- 
erator L is the adjoint of L in canonical correlations and 
the adjoint of the operator appearing on the right-hand 
side of the quantum master equation @ under the plain 
trace operation. For symmetric matrices 7,-fc, the only 
difference between the operators L and L is in the sign 
of the reversible term. 

The adjointness property ([12")) establishes a relation- 
ship or consistency between time evolution and canoni- 
cal equilibrium correlations. It may hence be regarded 
as a detailed-balance condition that ensures the proper 
symmetry of two-time canonical correlations. 



D. Detailed balance 

The usual construction of multi-time correlations re- 
lies on the possibility of introducing a Heisenberg pic- 
ture, that is, on the use of time-dependent operators av- 
eraged with a time-independent density matrix For 
a nonlinear master equation governing the evolution of 
the density matrix in the Schrodinger picture, the pas- 
sage to a Heisenberg picture is no longer possible. We 
are hence faced with a serious problem when we wish 
to study multi-time correlations. Even for the linearized 
master equation, one needs to be careful with employing 
the Heisenberg-like picture to define multi-time correla- 
tions. To discuss this problem in more detail, we start 
from the identity 

ti([A,B]p)=tr{(JC p A)[B,lnp}}, (14) 

which follows from Eq. ^ . When applied to the equilib- 
rium density matrix (J3]), we obtain the relation 

([A, B]) = i/3 (A; L rcv B) = (L rcv A; B) , (15) 

between equilibrium averages and canonical correlations, 
where £ rev and L rcv ~ —L rcv are the reversible contribu- 
tions to the super-operators L and L defined in Eqs. (jTOJ) 
and (fl3|) . After replacing A by e ct A, we obtain the rig- 
orous identity 

([e £t A, B}) = ij3 (l^ 1 A- b) . (16) 

For purely reversible dynamics, L — £ rov , we finally ob- 
tain 

([e^,B])=z/?|( e £t A;i?), (17) 

which is known as a fluctuation-dissipation relation be- 
tween a response function and a correlation function (see, 
for example, [1] or Eq. (4.2.18) of [IH). As it has been 
derived for from the fundamental Hamiltonian equation 
of motion, where A(t) — e ct A is the time-evolving opera- 
tor of the Heisenberg picture, the fluctuation-dissipation 
relation (fT7|) should be respected by all coarse-grained 
evolution equations, for example, the quantum master 
equation. 

For L L lcv , it is not allowed to define various 
two-time correlations naively by inserting e ct A because 
Eq. (|16[) would then imply a violation of the fluctuation- 
dissipation relation. This subtlety is known as the failure 
of the quantum regression hypothesis @, [ID, [H, 03 • It 
is related to the fact that an adjointness property like 
in Eq. (TT2l does not generally exist for the average of 
commutators, 

([LA, B})- ([A, LB}) = 

(L rev LA; B) - (A; L rev LB) ) = 

i/3 (A; (LL rcv - L rcv L)B) , (18) 
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which follows from Eqs. (|T5|) and (|T2|) . Adjointness, or 
detailed balance, in commutator averages can be guaran- 
teed only for purely reversible dynamics. 

E. Perturbation theory 

In the following, we are interested in frequency- 
dependent correlations rather than in the time-dependent 
correlation function (|17[) . We hence introduce the super- 
operator 

/>oo 

TZ(u>) = / e {£ - lu)t dt. (19) 
Jo 

If necessary for convergence, the frequency u> can have a 
small negative imaginary part. For nontrivially interact- 
ing systems, we would like to calculate such correlation 
functions by means of perturbation theory. We hence as- 
sume that the total Hamiltonian can be written as a sum, 
H = ij' ' + where for the super-operator £^ ob- 

tained for H( ' instead of H in the definition ([13")) (note 
that also K. involves H through the canonical equilibrium 
density matrix), but with the same coupling operators 
Qk, the frequency-dependent super-operator characteriz- 
ing the free evolution, 

/>oo 

U^(cj)= / e^-^dt, (20) 
Jo 

can be evaluated explicitly. The operator (w) is de- 
fined in the same way in terms of £^ \ These definitions 
imply the useful identities 

C<® + [T^H]- 1 = m + ^(u)]- 1 = iw. (21) 

We now consider the straightforward formal second- 
order perturbation series for the frequency-dependent 
canonical correlation, 

(U(u)A;B) = (k^\uj)A-B^ 

+ (n i0 HLu)C"JZ^(io)A;B) 

+ (tZ^\uj)£'TZ {0 \u)£'TZ W (u)A ; b) 
+ ... (22) 

where we have introduced the following pair of super- 
operators describing interaction effects, 

£' = £-£ {0 \ £'=£-£ {0 \ (23) 

Whereas the adjointness property (TT21 implies that the 
exact result is equal to (A; 1Z{uj)B), such detailed bal- 
ance properties do not exist for the individual terms of 
the perturbation expansion. The construction of a proper 
perturbation theory is even more subtle than the imple- 
mentation of the fluctuation-dissipation relation because 



the perturbation expansion involves multi-time correla- 
tions. In order to arrive at a more symmetric formu- 
lation, we need to find a suitable adjointness property 
involving £' and £' . From Eqs. (fT2"j) and (|2ip . we obtain 
the adjointness property 

U'n^(Lj)A;TZ^(u})B) + (k<®(w)A;B) = 
<^ (0) (cj)A; C'llW + (A; ft (0) (u)b) . (24) 

Moving £' to the other side is linked with moving TZ^> (w) 
to the other side in a lower-order term. The structure 
of terms of successive orders in perturbation theory be- 
comes coupled. After moving the first factor 1Z^ '(oj) in 
the highest-order term of Eq. (|2"2"j) , the adjointness prop- 
erty (|24[) allows us to rewrite the perturbation expansion 
in the more symmetric form 

(n(u)A;B) = (K {0 \uj)A;B) 

+ (£'K i0) (iu)A;K (Q HLu)B^ 

+ (k( Hlo)£'K<> \lo)A;£'kW{u)b) 
+ ... (25) 

The occurrence of super-operators acting on A and B is 
now balanced. 

In practical calculations, averages of commutators are 
more convenient than canonical correlations. By analogy 
with Eq. (|25l) we hence introduce the following correla- 
tion function guided by symmetry in time, 

C ab {lo) = ([1l^{w)A,B]) 

+ ([£"R {0) {u)A,1l (a \Lo)B]^ 

+ ( [K (Q) {u)£'n [Q) (u)A, £'K (0) (cj)B]j 
+ ... (26) 

For Hamiltonian dynamics, with £ = £ rcv , this pertur- 
bation expansion can be obtained in exactly the same 
way as Eq. ([231) because Eq. (TT51) provides the required 
adjointness property. For general master equations, with 
£ =/= £ rcv , we proceed as for the fluctuation-dissipation 
relation and postulate that the balanced perturbation ex- 
pansion ([26| still must be applicable. The validity of this 
postulate can be verified in the applications of Eq. (|2"6")l 
in Sees. |IVH and EEE 

A more rigorous approach to time-reversal and detailed 
balance properties in perturbation theory will be of par- 
ticular importance if one is interested in multi-time cor- 
relations, such as Wilson loops. We need a generalization 
of the discussion of time-reversal symmetry and detailed 
balance for the fluctuation-dissipation relation in Sees. 4 
and 6 of Q. For our purposes, however, the balanced 
second-order perturbation expansion (|26[) turns out to 
be fully satisfactory. 
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The averages in Eq. (|26|) are still for the interacting 
theory. In the following, we need the first-order expan- 
sion (see, for example, Eq. (4.1.10) of 15], or Eq. (5.1.9) 
of0) 

(A) = (A) (a) - ft ((A; F (1) )) (0) + . . . , (27) 

where the notation ((•; •)) implies that, in the evaluation 
of the free canonical correlation by means of Wick's the- 
orem, at least one contraction between A and ffW needs 
to be present. In other words, the contributions from a 
product of averages is suppressed so that we might call 
((■; •))'°) a canonical covariance. 

III. RENORMALIZATION GROUP 

Renormalization is sometimes perceived as a tricky 
toolbox to remove annoying divergences from quantum 
field theory. The present section emphasizes that the 
renormalization group should rather be considered as a 
profound tool to refine perturbation expansions which 
would actually be useless without renormalization. We 
first try to develop some intuition in the context of poly- 
mer physics where no divergences are present. We then 
present a few equations that provide practical recipes for 
refining perturbation expansions and allow us to calcu- 
late the famous ft function for the running coupling con- 
stants of renormalization-group theory. 

A. Intuitive example 

Plain perturbation theory is clearly not appropriate 
for problems involving a large number of interactions. 
However, if a problem exhibits self-similarity on different 
length scales, perturbation theory can be refined to ob- 
tain useful results by successively accounting for infinitely 
many interactions. This refinement may be considered as 
a kind of generalized exponentiation procedure guided by 
a renormalization-group analysis. 

Nice illuminating examples of refined perturbation ex- 
pansions can be found in the theory of linear polymer 
molecules. The beauty of polymer physics actually stems 
from the self-similarity of polymers If we model 

polymer molecules in dilute solution as linear chains of 
beads connected by springs, hydrodynamic interactions 
between the beads arise because each bead perturbs the 
solvent flow around it and, after propagation of the per- 
turbation, it affects the motion of the other beads. The 
bead friction coefficient determines the strength of such 
hydrodynamic-interaction effects. 

The beads of such mechanical polymer models, how- 
ever, are fictitious objects consisting of many monomers. 
If one uses larger beads, hydrodynamic interactions be- 
tween the smaller beads inside a larger bead have to be 
incorporated into the effective friction coefficient of the 
larger beads. This consideration offers the possibility of 



incorporating more and more interactions by passing to 
successively larger beads. Small increments in bead size 
implying only few interactions inside the larger beads can 
be handled by perturbation theory and renormalization- 
group theory allows us to accumulate a very large num- 
ber of interactions via many small steps. A number of 
static and dynamic properties of polymers in solution 
have been computed by these ideas [19h24| . The transfor- 
mation behavior between effective interaction strengths 
on different length scales contains important informa- 
tion, in particular, about the critical exponents associ- 
ated with self-similarity and the limiting value of the ef- 
fective coupling strength on large scales. With the large 
scale model, one can finally perform a perturbative cal- 
culation of any quantity of interest. A comparison with 
the scaling expressions for these quantities following from 
self-similarity allows us to refine the results of plain per- 
turbation theory even further (see, for example, the be- 
havior of various material functions at high shear rates 
discussed in [25|). 

Note that divergences are not an issue in the above dis- 
cussion of hydrodynamic interactions in dilute polymer 
solutions. They would only arise if, to establish contact 
with field theory, we considered the limit of an infinitely 
large number of infinitesimally small beads. In the next 
section, we show how the intuitive ideas of this section 
can be translated into tractable equations. 



B. Basic equations 

Let us introduce a small length scale i, which could be 
a bead size, a lattice spacing, an inverse cutoff for mo- 
mentum or frequency, or the characteristic length scale 
of dissipative smoothing. If some model of interest con- 
tains a dimensional coupling constant A, say the strength 
of some interaction, the proper choice of which typically 
depends on £, we first introduce the dimensionless cou- 
pling constant 

\{i)=l*\{l), (28) 

with suitable exponent e obtained from dimensional anal- 
ysis. The intuitive ideas described in the preceding sub- 
section are implemented by constructing a perturbation 
theory of the rate of change of A(^) with £ rather than for 
A(£) itself. We hence assume that there actually exists 
a perturbation expansion of the function describing the 
rate of change of the dimensionless coupling constant, 

P$) = -l% (29) 

that is, for the standard ft function for the running cou- 
pling constant of quantum field theory. We always dis- 
play the ft function with its argument to avoid confu- 
sion with the inverse temperature. If the free theory re- 
mains free on all length scales, ft(0) = 0, the most general 



6 



second-order expansion of the /3 function is given by 



/3(A) = -aX 1 - — , (30) 



where the parameters a and A* remain to be determined. 
According to Eq. (f2U| , the second-order perturbation the- 
ory (|30|) for /3 implies the following nontrivial dependence 
of the running coupling constant on the length scale, 



X(£) 



A* t° 



(a 



Note that the combination 



A* 
c 



£ e X(£) 



X* 



(31) 



(32) 



is a constant. Any invariant property P must be some 
function of this constant and, as a polynomial in A*/c 
and eventually in X(£), must possess a second-order per- 
turbation expansion of the form 



Po + Pi 



XI 



e e x 
~x* 



p f( a -e) ^ 



(33) 



We will later find the identity a = e so that the per- 
turbation expansion (|33[) takes an even simpler form. 
Then, Eq. (|3"Tj) shows that A(0) corresponds to A*/c, and 
X{£) = £- e X* for large £. 

Of course, these arguments can be generalized to con- 
struct the most general form of higher-order perturba- 
tion expansions of observable quantities. Generalizing 
Eq. (f3"2")> . a polynomial expansion of the /3 function leads 
to a nonpolynomial form of the constant, 



A = Cr'expi-e 



A dX' 



/3(A') 



(34) 



where the integration constant C is to be chosen such 
that A = A for small £. The integral of the rational func- 
tion 1//?(A') can actually be performed in closed form, 
but only an expansion in A is required. Any invariant 
quantity P must be a function of A and hence possesses 
a perturbation expansion of the form 

P = P + PiA + P 2 A 2 + P 3 A 3 ..., (35) 

where the expanded form of Eq. (|34l) has to be inserted, 

A = A + B 2 £ f X 2 + B 3 £ 2e X 3 .... (36) 

The combined expansion is of the form 



P = P + PiA + (P 2 + PiB 2 £ e )X 2 
+ (P 3 + 2P 2 B 2 t + PiB^X 3 



(37) 



Note that this expansion contains the various corrections 
from the small length scale £. The coefficients Pj can 



easily be read off from the terms of order £° of such an 
expansion. The coefficients Bj characterizing the /? func- 
tion can then be read off in various places; in the language 
of Feynman diagrams, this possibility corresponds to the 
occurrence of the same subdiagrams in infinitely many, 
increasingly complicated diagrams. The consistency of 
the results and their independence of the particular quan- 
tity P expresses the renormalizability of the theory. 

Note that the dependence of A on £ is obtained from 
the analysis of the perturbative prediction of large scale 
properties. The parameters (X(£i),£i) and (X(£ 2 ),£ 2 ) 
with £ 2 > £\ imply the same large scale properties (on 
length scales large compared to £ 2 ), but the model with 
parameters (X(£ 2 ),£ 2 ) typically cannot accommodate all 
properties of the model with parameters (A(£i),^i) on 
length scales of order £ 2 . In other words, the underlying 
model serves as a minimal model for a universality class 
without being close to some fixed point model. If we 
want to refine our perturbation theory for some observ- 
able P, however, we only have the option to translate the 
coupling constant A from small scales to some physical 
length scale, say an inverse mass m, and to use the cor- 
responding translated A as A in the expansion (|35p . An 
even simpler refinement is obtained if we assume that, 
at the length scales of interest, A has reached its fixed 
point value A* with (3(X*) = 0, which implies the choice 
A = m e X* in the expansion (l3"5"j) . or A = m e X* and £ — 
in the expansion (|37p . 

Keeping the large scale physics invariant when letting 
the length scale £ of the underlying model go to zero im- 
plies that the underlying model has to approach a crit- 
ical point with diverging correlation length in units of 
£. This well-known relationship between field theory and 
critical phenomena, including the calculation of critical 
exponents, has been discussed extensively, for example, 
in the nice classical review article [26J . 



IV. 



With the tools of quantum dissipation at hand, we 
can now formulate quantum field theory with small-scale 
smoothing by a dynamic friction mechanism. Our ap- 
proach may be considered as a generalization of the 
canonical quantization procedure to include dissipation. 
We present the free field theory and construct the prop- 
agator and a four-point correlation of the interacting 
theory in second-order perturbation theory. In a fi- 
nal step, we refine perturbation expansions by means of 
renormalization-group theory and construct the /3 func- 
tion. 



A. Quantization procedure 

Quantization procedures are traditionally based on a 
canonical Hamiltonian formulation of the evolution equa- 
tions for the underlying classical systems. In the canon- 
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ical approach to quantum mechanics, rooted in Dirac's 
pioneering work, the canonical Poisson brackets of classi- 
cal mechanics are replaced by the commutators of quan- 
tum mechanics [27j |. This procedure has been adapted 
to quantum field theory (see, for example, Sees. 11.2 and 
11.3 of [I| or Sec. 1.8 of ;2<|). Even in the path-integral 
approach to quantum field theory, the justification of the 
proper action needs to be supported by the canonical 
approach (see, for example, the introduction to Sec. 9 of 
|30|). Understanding its Hamiltonian structure is hence 
crucial for quantizing a system. 

In classical mechanics, on has the choice between the 
Lagrangian and Hamiltonian formulations; both formu- 
lations can also be used for classical field theories (see, 
for example, Chaps. 2, 8, and 12 of [3l[). The Lagrangian 
approach is based directly on the variational principle for 
the action which is obtained as the time integral of the 
Lagrangian. In contrast, the equivalent Hamiltonian ap- 
proach needs two structural elements, the Hamiltonian 
and a Poisson bracket required to turn the gradient of 
the Hamiltonian into the vector describing time evolu- 
tion. The canonical Poisson structure, which in the non- 
degenerate case is also know as symplectic structure, is 
the key to formulating the proper commutators in the 
quantization procedure. It is important to note that the 
Hamiltonian and the Poisson structure are two separate 
elements. In particular, once one has identified a fun- 
damental Poisson structure, one can easily change the 
Hamiltonian. To understand the equivalence of the La- 
grangian and Hamiltonian formulations of classical field 
theory, one should realize that the existence of a non- 
degenerate Poisson structure is crucial to establish an 
underlying variational principle [32| . 

If one changes the system, one needs to change the La- 
grangian and hence one needs to begin from scratch in 
the Lagrangian approach. In the Hamiltonian approach, 
one needs to change the Hamiltonian but one may keep 
the Poisson structure, where the latter is the key to quan- 
tization. Of course, one might ask whether we actually 
want to change the system. This indeed happens when 
we consider the noninteracting and interacting systems 
in parallel, say for constructing a perturbation theory. 
The logical separation of canonical Poisson structures 
and Hamiltonians is also very useful when one wishes 
to quantize systems involving nontrivial constraints. In- 
deed, the proper handling of the constraints resulting 
from gauge invariance is a major obstacle to quantizing 
the Yang-Mills and gravitational fields. Note that the La- 
grangian and Hamiltonian approaches are closely related 
to path-integral and canonical quantization, respectively. 

As noncquilibrium thermodynamics is built on a 
Hamiltonian formulation of reversible dynamics, the 
quantization of dissipative systems is obtained as a gener- 
alization of the canonical quantization procedure. As an 
additional step, one needs to formulate the proper form 
of the frictional coupling to a heat bath in the quantum 
master equation. 



B. Fields and Hamiltonian 

We consider a scalar quantum field (p(x) and its canon- 
ical conjugate tt(x) in d space dimensions. Throughout 
this paper, we denote the dimensions of space and space- 
time by d and D = d + 1, respectively. The bosonic 
field operators <p{x) at all positions x commute among 
each other, and so do the conjugate operators tt(x). The 
only nontrivial commutation relations are of the canoni- 
cal form 

[(p(x),ir(x')] =iS(x-x'). (38) 

We assume that the total Hamiltonian of our system is 
given by 

H = ||i[7r 2 + (V^) 2 + mV]+-^A/}^, (39) 

where the quadratic contribution in tt(x) and f(x) de- 
scribes a free massive scalar field with mass parameter 
m and the quartic contribution is the interaction term 
of the ip 4 theory. The parameter A, often referred to as 
a coupling constant, describes the strength of the inter- 
action. A dimensional analysis of Eq. (|3T)f shows that 
ip 2 m 1 ~ d and Xm d ~ 3 are dimensionless. 

The self-adjoint position-dependent canonical field 
variables can be expressed in Fourier representations of 
the form (see, for example, Sec. 1.8 of 29] or Sec. 12.1 of 
.28]) 




The only requirement for the function u>k occurring in 
this representation is that it assumes equal values for k 
and — k. With these Fourier representations, we have in- 
troduced the adjoint operators a' h and a k creating and 
annihilating field quanta of momentum k G Mr, respec- 
tively, as primary variables. All creation operators com- 
mute among each other, and so do the annihilation oper- 
ators. The only nontrivial commutation relations for the 
boson creation and annihilation operators are 

[a k ,a{,}^S(k-k'). (42) 

We assume that the collection of all states created by 
multiple application of all the operators a/ k for all k £ K d 
on a ground state (which is annihilated by any a k ) is 
complete. The full Hilbert space factorizes into spaces 
obtained by repeated application o f at for each mode k 
(see, for example, Sees. 1 and 2 of [33[ or Sees. 12.1 and 
12.2 of [28] for more details on the construction of such 
Fock spaces). The field quantization based on Eq. (|4"2")l 
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for creation and annihilation operators is an equivalent 
alternative to the canonical quantization procedure based 
on Eq. (|38|l. 

We next write the total Hamiltonian as the sum H = 
fj(o) + fji 1 ) anc j express the quadratic free Hamiltonian 
and the quartic interaction term H^' in terms of 
creation and annihilation operators. Neglecting an irrel- 
evant constant contribution, the free Hamiltonian can be 
expressed in the simple form 



H 



(o) _ 



to k a k a k d d k, 



(43) 



provided that the momentum-dependent frequencies w k 
are given by the relativistic dispersion relation 



«2 



k 2 



(44) 



Once more neglecting constant terms, the interaction 
term of the if 4 theory is given by 



A 1 



--- rl d h 



+ k 2 + k 3 + ha) 



96 (2w) d 

(a- kl a-k 2 a-k 3 a- ki + 4aj Ci a_fc 2 a_fc 3 a_ fc4 
+ 6 4 1 a L a -fc3 a -fc4 

+ 4 4 1 «L a L a -fe4 + 4i4 2 a L a L) 

, f d d k ( n t t t \ 

+ A2 / --- [a k a- k + 2a k a k + a k a*_ k j , 



(45) 



where we have chosen to use the normal ordered form of 
that is, all creation operators are moved to the left, 
all annihilation operators are moved to the right. Nor- 
mal ordering will be very convenient for the subsequent 
calculations. By inserting the Fourier transform (|40[) into 
the Hamiltonian (|39p , we would actually obtain = 



1 1 



4 (27r) d 



d d q 
2lj„ 



(46) 



where the integral I\ and its Lorentz invariance are dis- 
cussed in Appendix lAl [see Eq. (|A1[) ]. 

We have introduced the generalization from u/fc to oj k 
with the idea to modify the frequency at very large |fe| 
to guarantee the convergence of all integrals occurring in 
intermediate steps of the calculation and to allow us a di- 
rect comparison with more traditional regularization pro- 
cedures. The details of this regularization at large wave 
vectors are irrelevant to the present approach. Note that 
the passage from u>k to Q k does not change the Hamil- 
tonian structure (only the Hamiltonian), but it destroys 
the Lorentz invariance of the system. 

It is very convenient to allow for values of z that are 
different from the value occurring naturally from the 
normal-ordering procedure, z no given in Eq. (|46[) . To fix 
the value of z in a more convenient way, we look at the 



most fundamental correlation function (|20[) obtained for 
A = a k and B = at,, that is, for a field quantum created 
with momentum k' and later annihilated with momen- 
tum k. We introduce the correlation function C k {u> 2 ) by 



c 



» 



2u 



C k (u 2 )5{k-k'), (47) 



where momentum conservation has been taken into ac- 
count and, in view of its symmetry in u, a dependence of 
C k on uj 2 has been anticipated. By choosing z such that 
the moment condition, 



C (w 2 



dC (u 2 ) 



j 2 =0 



duo 2 



(48) 



j 2 =0 



is fulfilled, we have a well-defined relationship between 
the mass parameter m and a physical correlation func- 
tion of the interacting system in the limit of small fre- 
quencies and zero wave vector. With this choice of z, 
the physical meaning of mass is not affected by inter- 
actions. This leads to a nicely organized perturbation 
theory in A that can be refined very conveniently by a 
renormalization-group analysis and avoids the usual need 
for mass renormalization. Allowing for z ^ z no results in 
what is usually achieved by a counterterm. We assume 
that z possesses a perturbation expansion in A. It turns 
out below that this perturbation expansion of z can be 
constructed very easily, in particular, if we rewrite it in 
the equivalent form 



Co(Q) 



______ 

n 2 duj 2 



{lo 2 - m 2 ) 2 C {w 2 ) 



(49) 



We further introduce the factor Z by the condition 

C o (0) = ZC£ 0) (0), (50) 

so that we can preserve the normalization in the presence 
of interactions by rescaling the field operators or, equiva- 
lently, the creation and annihilation operators. Equation 
(|50| offers an explicit recipe for calculating Z. In partic- 
ular, a perturbation theory for the two-time correlations 
on the left-hand side of Eq. ((50)) directly provides an ex- 
plicit perturbation expansion for Z. Note that, in the 
thermodynamic approach, a change in the normalization 
of the fields in the Hamiltonian can trivially be absorbed 
in a change of the coupling constant A and an overall 
factor to be included into the definition of temperature. 

C. Friction mechanism 

We now apply the linearized thermodynamically con- 
sistent master equation (|10p to the treatment of quantum 
field theory. With the coupling operators 



Qj, Qk ^ vVO), vV(y), 



(51) 



which favor dissipation and smoothing of small-scale fea- 
tures through the occurrence of the Laplacian V 2 as the 
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simplest scalar differential operator, and the local relax- 
ation rates 



j jk H> 2-/6{x - y), 



(52) 



we arrive at the preliminary quantum master equation 



d P ■ r- rn / ,d, Ik 

x /C _1 [afc t + a_fc,/C [a_fc' + afc,/5]], 



(53) 



where 7^ — 7|fc| 4 with the friction parameter 7. Note 
that the occurrence of |fc| 4 is natural because the scalar 
coupling operator, which involves the Laplacian, occurs 
quadratically in the quantum master equation. 

The quantum master equation (|53[) is obtained un- 
der the assumption that only the Laplacian of the fields 
should be used for the dissipative coupling, not the Lapla- 
cian of their canonical conjugates. A more symmetric 
coupling to the fields and their conjugates leads to the 
simpler final equation 

x r^KUh.pE + K^K^]]). (54) 

Such a symmetric coupling also arises for the damp- 
ing of electromagnetic field modes inside a cavity (see 
Eq. (3.307) of [3]). For single creation or annihilation 
operators, Eqs. @, © and (|4"2")l lead to the remarkably 
simple results 



and 



Cal 



.Tfc 



[a k ,H] 



(55) 



(56) 



The quantum master equation (|54|) is our fundamen- 
tal equation of quantum field theory obtained by a gen- 
eralization of the canonical quantization procedure for 
dissipative systems. In addition to the usual reversible 
evolution given by the commutator with the Hamilto- 
nian, local degrees of freedom are damped by an irre- 
versible friction mechanism expressed in terms of a dou- 
ble commutator. The cutoff of a quantum field theory 
is realized as a spatial smoothing achieved by a dissipa- 
tive dynamical mechanism where smoothing over short 
distances takes place very quickly. 



D. Free field theory 

In this section, we consider the calculation of equilib- 
rium averages and the evolution equations for the free 
theory with A = 0. To generalize Wick's theorem for 



evaluating complicated moments to the case of finite tem- 
peratures, we use the ideas of the chapter "field theory 
at finite temperature" (Chap. 7) of [33[ (see also the the- 
ory of temperature Green's functions in Sec. 5.6 of [Hj]). 
Note that equilibrium averages are not affected by the 
fact that we use a quantum master equation with dissi- 
pation rather than a purely reversible Hamiltonian evo- 
lution. 

In the spirit of Eqs. (24.32) and (24.36) of jH, we 
obtain for averages performed with the free Hamiltonian 
for an arbitrary observable A: 



(0) 



akAy = — — 



ak,A] 



AO) 



(57) 



Wick's theorem (f57| is the working horse for evaluating 
the free averages containing an increasing number of cre- 
ation/annihilation operators and hence plays a crucial 
role in perturbation theory. 

In the following, it is often convenient to consider ob- 
servables of the normal ordered form, 



■ a k i «fci • ■ • a kK ■ 



(58) 



In particular, we are interested in the free evolution of 
such a normal ordered operator A. We study the ingre- 
dients to the reversible and irreversible contributions in 
the evolution equation separately. 

The reversible part of the free evolution is obtained 
from the commutation relations 



: uj k a k . 



Ofcl 



-uJka k - 



(59) 



These commutation relations imply the more general re- 
sult 



[H<-°\A]=u> A A, 



where the frequency 



K 



3=1 3=1 



(60) 



(61) 



has been associated with the operator A. 

By following the change with it, we obtain as a useful 
consequence of Eq. (fSU)) 



A, 



(62) 



for arbitrary complex u. The normal ordering of the op- 
erators in Eq. (|5"5|) is actually irrelevant for the above 
results (fFJO")) and (|6"2"|) . We further obtain 



(63) 



K,^ A = w{uj a ) Ap 1 - ^ = w(—oja) 



where p^ is the equilibrium density matrix ^ for the 
free field theory and the function w(tu) is defined as 



(64) 
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Equation (|63|) follows from Eqs. (|2|) and (|62[) by explicit 
integration over u. From Eq. (|63p we further obtain 

(A; B) {0) = «;(-wa) (AB) (0) = w(oj b ) (AB) {0} , (65) 

and, with the additional help of Eq. (l62l) , 

([A, B]) {0} = (A; S) (0) = 0u> B (A; B) W , (66) 

so that all kinds of free equilibrium averages are related 
by simple factors. Note that Eq. (|66|) can also be ob- 
tained as a simple special case of the much more general 
result (|15p for the interacting system. 

The irreversible part of the evolution super-operator 
of the free theory can now be evaluated by means of 
the identity (f6"5|) . As a first step, we obtain the useful 
identities 

k.w\a,k^b] = 

Puja ^a)w(u Jb ) ba+ W ( UB ) 

W(UJ A +UJ B ) w{<jJA+Ub) 

w(uj a )w(uj b ) w{uj B )e-^ A 

PUA-, ; V AB + — ■ r[A,B\. (67 

w{uj A +u B ) w{uj a +uj b ) 

After using Eqs. (g>0) and (|67j) in Eq. ([gljl. the total free 
evolution super-operator resulting as the sum of the re- 
versible and irreversible contributions is found to be 



L {0) A = —iloaA — / d d k 



Ik 



w(uj a ) 

x I w(io A - LJ k )w(u) k )a k ^[a k ,A] 



+ w(u A + oj k ) w(-u k ) [A, a k 1 ]a k 



-w(uj a - uj k ) [[a fc , A],a fc t ] 



- -^-w(uj a + uj k )[a k ,[A,a k f ]]j. (68) 

Note that the double commutators in the last two terms 
in Eq. ([55]) are actually equal. According to the first two 
terms under the integral in Eq. (|68[) . the decay rate i A 
for an operator A is given by 



7A 



y~l ik'W(uj k <. , uj A - w fe / ) 



K 



^lk j W{-w kj ,uj A + ^ kj ), (69) 



with 



W(u,u') = 



(70) 



w(u + UJ 1 ) 

The function W has the useful symmetry properties, 
W{w,u') = W{w',uS), W{uj,uj') = W{-u,-ui'). (71) 



By simplifying the result (1681) for normal-ordered prod- 
ucts (|58[) of creation and annihilation operators, the fun- 
damental evolution operator of the free theory can now 
be rewritten in the compact form 



£^A = -(iLo A + lA )A + r A , 



(72) 



where the term T A consists of all the contributions ob- 
tained by deleting one creation and one annihilation op- 
erator from A. The explicit normal ordered form of F A 
is given by 



r A = 



d d k-^ -\w(ui A + uj k ) 

f3uj k w(uj a ) I 

+ e-P Uk w(lo a - w fe )} [a k , [A, a],] 



(73) 



When exponentially small terms are neglected in the low- 
temperature limit, Eq. ([73"]) for u) A ^ can be simplified 
to 



Ik \uj A \ 



f3uj k (ujk + \u>a\) 



[a k , [A, at] 



(74) 



Some examples of T A for normal ordered products of up 
to three creation/annihilation operators are given in Ap- 
pendix [5] 

From the identity (|2"Tj) and the explicit expression ([72"]) , 
we obtain the result 



K(°\uj)A 



1 



ius + iuj a + 7^ 



A + n {0) (uj)r A , (75) 



by which the calculation of TZ}- \<jj)A can be reduced to 
successively simpler products of creation and annihilation 
operators. We similarly have 



n w (tu)A 



i 



ius — itu A + j A 



A + 7^ )(w)rJ . (76) 



E. Propagator 



We are now ready to analyze the fundamental correla- 
tion function C k (uj 2 ) introduced in Eq. (|47p for an inter- 
acting system. If we introduce the equilibrium averages 



{{a k ,H(%al])=X k 5(k-k'), 



(77) 



and 



[RF>(u)[a h ,HW], [HW, aU]) = i Y k {u) S(k - k'), 



(78) 

and use the results ([55]). ([56) . (1751) . and ([751), tnen 
Eqs. fl25J and (gTJ lead to 

C k (uj 2 ) 



+ 



oj 2 - (u> k - ij k ) 2 
2i(uj k - ij k ) 2 
ul[u 2 - (uj k ~ ilk) 2 ] 2 



u k X k + (u k - i"f k )Y^(ui) 



i{uj k - ilk) 2 [u 2 + (wk ~ ilk) 2 ] Y k (cj) 
ojI[uj 2 - (uJk - ilk) 2 ] 2 w 



(79) 
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FIG. 1: Feynman diagrams potentially contributing to the 
propagator in second-order perturbation theory. 



where Y k + (oj) and Y k ~ (oj) are the symmetric and anti- 
symmetric contributions to Yk(oj). 

The exact results for X k and Yk(oj) for any tempera- 
ture are given in Appendix [DJ Equations ([Dip and (|D3|) 
allow us to see explicitly how the limits of low temper- 
ature and weak friction can be performed. Neglecting 
exponentially small terms at low temperatures, we find 



X 



fc — — , 



(80) 



and 



Y h (u) = 



r,2 



A 2 1 f d d k 1 d d k 2 d d k 3 



96w fe (2tt) 



2d 



with the damping rate 



<5(fei + k 2 + k 3 + k) 
(81) 



P 

lk 2 



Ik, 



i 7fc 3 

oj k2 (Wfcj + Wfc 3 ) w fe;3 (w fcl + uj k2 ) 

and the rational kernel function 

1 1 



OJ + OJ — l"f OJ — OJ — IJ 



(82) 



(83) 



The occurrence of 17 in Eq. (|83[) implies specific rules 
for the treatment of the poles in the integrations of 
Eq. (f5Tj) . The 17 hence plays a similar role as the ie 
in the usual approach (see, for example, Eqs. (|A1|) . (|A8[) 
and the comments offered there). The usual ie resolves 
infrared problems by clarifying causality issues. In our 
approach based on irreversible equations, there is a nat- 
ural arrow of time and a large scale decay of correlations 
that leads to convergent integrals. Therefore, the small 
friction coefficient simultaneously resolves infrared prob- 
lems by ensuring a decay of correlations on large scales 
and ultraviolet problems by dissipative smoothing of the 
fields on small scales. 

As the averages in Eqs. (1771) and ([781 need to be evalu- 
ated by means of Wick's theorem, we can introduce Feyn- 
man diagrams to keep track of the contraction structure 



of the various contributions. The contribution X k cor- 
responds to the diagram in Fig. [TJa), whereas the con- 
tribution Yk(w) consists of two terms represented by the 
Feynman diagrams in Figs.fljb) and (c). The diagram in 
Fig. QJa) looks like a first-order contribution, as also the 
factor of A in Eq. (f8U[) suggests. Note, however, that this 
contribution stems from the quadratic part of the inter- 
action Hamiltonian (|45[) arising in the normal-ordering 
procedure and that it hence contains a factor of z, which 
is shown as an important reminder in the Feynman dia- 
grams. This factor z still needs to be determined from 
Eq. (|4"8"j) or (14*§1) . It turns out that z itself is of order A so 
that the diagrams in Figs.QJa) and (b) actually represent 
second and fourth order contributions, respectively. The 
contribution in Fig.[ljc) clearly is a genuine second-order 
term as suggested by the factor A 2 in front of the triple 
momentum integral in Eq. f|81|) . 

In the limit of vanishing friction, we can set 7fe = in 
Eqs. (|7J) and ([8"T]). Equation {7JJ becomes 



{oj 2 - oj 2 k f C k {oj 2 ) = i{oj 2 -oj 2 k )+2iu k X k + Y+(oj) 



(84) 



In simplifying Eq. ([8ip we need to be more careful with 
7fc 1 fe 2 fe 3 because this quantity keeps all integrals finite, 
which expresses the regularizing character of the smooth- 
ing friction mechanism. To simplify the kernel function 
occurring in Eq. (|8ip. we consider the following expan- 
sion in to, 



OO 

E 

n=0 



1 



(u> + ij) n+1 (oj — 17) 



ra+l 



(85) 



After neglecting the regularizing friction mechanism in 
all convergent integrals (where, for simplicity, we assume 
d < 3, so that the integrals converge for n > 2), we obtain 
the simplified expansion 



1 



1 



17 oj 

,2n 



17 



(OJ - 



OJ 



il) 2 
(86) 



We next analyze the terms proportional to u in Eq. I[86p . 
Note from Eq. ([82]) that the regularizing parameter is 
7//3, and dimensional analysis hence implies that the cor- 
responding smoothing length scale is given by 



(87) 



where the factor 2 is introduced for later convenience. 
Each of the resulting two contributions to Y k propor- 
tional to oj diverges as £-( 2d ~ 5 ) (for d < 5/2, the terms 
are convergent and cancel each other in the limit of zero 
friction). For symmetry reasons, the divergent terms (for 
5/2 < d < 3) cancel each other and the corrections, when 
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expanded in terms of the dimensionless parameter £m, 
vanish in the limit of small friction. We can hence rewrite 
the kernel function (l86l) as 



r(u),ui,i) 



2u 



■T 



00 1 , 2n 
. 2 V_ 

ra=l 



(88) 



implying YjT(u)) = Yk(ui) and (u>) = for the sym- 
metric and antisymmetric contributions to Yfc(w). At the 
expense of losing the convergence of the integral of each 
individual term, Eq. (|88[) can formally be re-summed into 
a simple rational function of lu 2 and a constant term with 
prefactor j 2 , 



r(w,w,7) 



2Cu 



27 2 



lu(lu 2 + 7 2 ) 



(89) 



After simplifying the kernel function, we are now in a 
position to discuss the propagator in more detail. For the 
second-order expansion Ya(u) = Y+Y'u 2 , the conditions 
(|49|) and (|50|) for the correlation function (|84|) imply 



and 



X = —(Y + m 2 Y'), 



Z = 1 + 2mY'. 



(90) 



(91) 



Note that the value of z following from Eqs. (|9"0"J). ([5D]>. 
and (f8"T1) with the kernel function for which we ob- 
tain the straightforward result 



A 1 f d d k 1 d d k 2 d d k 3 8(ki + k 2 + fc 3 ) 
48 (27r) 2d J Q kl ti k2 tik 3 w fel + w fc2 + w/c 3 
(o; fcl +^ fc2 +(Jfc 3 ) 2 +™ 2 

(W fel +Wfe 2 +W fc3 ) 2 +7fc lfe2 fc3' 



(92) 



turns out to be small, more precisely, only of first order 
in A. In terms of the integrals introduced in [34j and 
reproduced in Eq. (IA7|) . we find 



and 



Z=l + i(Am d - 3 ) 2 / 3 , 



z = y 2 m 2d -\h-h) 



(93) 



(94) 



The result (|93j) for the normalization factor Z coincides 
with the one given in Eqs. (11.64) and (11.69) of [34l |. 
The quantity z is usually taken into account through 
mass renormalization. More precisely, we have verified 
explicitly that Z [1 + 2 A (z — z nn )/m 2 ] coincides with the 
mass renormalization factor Z m of [34J . 

For presenting the results of the perturbation theory, it 
is convenient to "amputate the external legs" associated 
with free propagators (see, for example, p. 55 of 
We hence introduce the usual amplitude 



(u, 2 



C k (uj 2 



(95) 



From Eqs. (|84|) and (|9Tj) . we obtain for the second-order 
perturbation series 

M k {u: 2 ) = -i(oj 2 -ujl)-2i\uj k Y k (u})~mY~(L} 2 -k 2 )Y' 

(96) 

with nicely Lorentz invariant correction terms due to the 
constants Y and Y' . The Lorentz invariance of 



WfcYfc(w) = - 



A 2 



1 



d d k 1 d d k 2 d d k 3 



96 (27r) 2d J o) fcl a) fe2 a) fe3 
x 5{ki +k 2 + k 3 + k) 

x r(ui,ui kl + oj k2 + uj k3 ,~f klk2k3 ), (97) 

is less obvious. Note that we have assumed u k = Lu k 
for the externally fixed momentum (regularization is only 
relevant in integrals). To establish the Lorentz invariance 
of WfcYfc(w), we need to pass from space to spacetime inte- 
grals. If we neglect all regularization effects by assuming 
also Cj k . — ui k . and by considering the limit of vanishing 
friction, 



2uj 



(98) 



A 2 /fc(w 2 )/12 in terms of the integral 
9j. Equivalent forms of the integral 



we find uj k Yk(uj 
defined in Eq. 

7fc(w 2 ) are given in Appendix [A"l The representation in 



Eq. (|A8j) demonstrates the invariance of uj k Yk(uj), clearly 
exhibits the structure of a product of three free space- 
time propagators, and reproduces the standard result for 
Mk{oo 2 ) as, for example, given in Eq. (III. 3. 2) of p9l| . 
Equation (|A6I) finally leads us to the Eucledian time ver- 
sion of the second-order perturbation series for the prop- 
agator given in Eqs. (11.36) and (11.38) of [H. 

Regularization mechanisms often destroy the Lorentz 
invariance, although a covariant formulation of friction 
within a not manifestly covariant thermodynamic frame- 
work may be possible (see, for example, [35- 37]). After 
neglecting the regularization mechanism for establishing 
Lorentz invariance, we next look at Eq. (IMj) or at Eq. (|9"?| 
with the kernel function ([88)) to see how integrals be- 
come regularized. In the traditional approach, one has 
7feifc 2 fc 3 — an d the convergence of integrals is achieved 
by an enhanced increase of the frequencies uj k at large 
wave vectors kj. In the present approach, we can use 
cu kj = uj kj because the occurrence of l klk2k3 keeps all 
integrals finite. 



F. Four-point correlation 

Our next goal is to obtain the j3 function for the cou- 
pling constant from a second-order perturbation theory 
of the form (|33[) . As a first-order term is missing in the 
perturbation series (|96l) for the propagator, we cannot 
determine the fixed-point value A* of the dimensionless 
coupling constant. We hence consider a four-point corre- 
lation which, to leading order, is given by the interaction 
strength A of (p 4 theory. 
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FIG. 2: Feynman diagrams contributing to the interaction 
vortex in second-order perturbation theory. 



As a second example, we hence apply our basic pertur- 



The 



bation formula (1261) to A = a kl a k2 and B 
result is of little direct interest because two particles are 
created at exactly the same time and later annihilated 
also at exactly the same time. Nevertheless, we will be 
able to produce a perturbation series from which we can 
read off the parameter A* and the j3 function (IBTfl) . It is 
actually sufficient to do the calculations in the limit of 
vanishing frequency u and wave vectors kj, k'-. 

According to Eq. (|T6|) , the zeroth-order contribution in 
Eq. (|26p factorizes into two propagators. After subtract- 
ing a suitable product of propagators, we are interested 
only in contributions corresponding to connected Feyn- 
man diagrams. As we restrict ourselves to the limits of 
vanishing frequency and friction (except for the regular- 
ization effect in integrals), we consider the property 



T = ([C' I0V A,B]) + (0 o Ho)L' Iev A,£' iev 



B] 



(99) 



With the explicit form of [[A, H<-%B] given in Eq. (jC9]) . 
we obtain 



([C' rev A,B]) = -i([[A,H^],B] 



-iX 



T 



with 



T = 



1 <5(fei + k 2 -k[- k' 2 ) 



(100) 



(101) 



where we have omitted all terms that vanish exponen- 
tially at low temperatures and all products of two prop- 
agators corresponding to particles that do not interact 
with each other. Note that no second-order correc- 
tions resulting from the expansion (|27p have survived in 
Eq. (fTUOl . 

According to Eqs. (|C3I) and (|C4[) . there are two types 
of contributions associated with the second-order term 
in Eq. (|59"|) . They can be represented by the Feynman 
diagrams shown in Fig.[5J For the contribution associated 
with Fig. [2ta), we have the explicit formula 



^ (0) (0) [a kl ,[a k2 ,H^]],[[H^\al { ],al, 2 
A 2 T f d d qi d d q 2 



32 (2ir) d 



■S(ki + k 2 + qi+ q 2 ) 



i(uj q 



•To. 



<2l " — Q2 

1 



iKi+w«,)-7 t iB t 



(102) 



where the expressions (IC5[) . (|C6[) and Wick's theorem 
(|57[) have been used to arrive at this integral expression. 
A typical contribution associated with Fig. [5fb) is given 
by 



ft (0) (0) 
A 2 T 



a kl ,H^]a k2 ,al [H™ 



d d qi d d q 2 



32 (2tt) c 



0J qi UJq 2 

8(qi + Q2 + ki - k[) 



J Q2 



+ Wfc^ +UJ k2 ) + 



; (103) 



where the evaluation is based on Eqs. (|C1|) and (|C2I) . 
In total, there are four contributions of the type (|103[) 
because fci and k 2 can be exchanged, as well as k[ and 
k' 2 . All other terms are exponentially small. 

Upon setting k\ — k 2 = k[ = k' 2 = in the integrals 
of Eqs. (fTU2")l and (1TU3]) and inserting them into Eq. ([55]). 
we obtain the simple perturbation expansion 



ImT 



T 



A + 



A 2 



1 



4 (2tt) 



d d q 



+ 



d d q 



(2ir) d J 2w»(w, + m)(l+7j' 2 ) 



with 



2 7 q 4 



f3uj q {uj q + 2m) 



(104) 



(105) 



In evaluating 7a- 51 a_, 2 , , and 7„_ gi0 _, 2 „ fc , Qfc2 , we 

have once more neglected exponentially small terms at 
low temperatures. In the definition of the four-point cor- 
relation, we have not included a factor 1/Z 2 because it 
would affect only third-order terms. 

For d < 3, the integrals in Eq. (|104j) are nicely con- 
vergent without any regularization. For u q = u q and 
7 = 0, these integrals have been expressed in terms of 
the r function in Eqs. (|A4I) and (|A5I) . By evaluating the 
leading-order corrections resulting from small values of 
the friction parameter 7, we arrive at 



ImT = - 

4 



A + Ja 2 m d - 3 (I{ +2J?) 



2 e 3-d 



(2tt) c 



qd d q 
1+q 4 



(106) 



where the smoothing length scale I has been defined in 
Eq. (|87p . We finally have a perturbation expansion with 
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FIG. 3: Fixed point value A* of the dimensionless coupling 
constant resulting from second-order perturbation theory as 
a function of space dimensionality d. 



a first-order term to compare to the general form (j3"3")l 
and we find a perfect match in structure. After using 
Eq. (IA2|) . the remaining one-dimensional integral can be 
evaluated in closed form (for example, with 3.241.2 of 
38] ) . From the comparison of the general form (|33p with 
our perturbation expansion (|106p . we then obtain the 
exponents 



a = e = 3 — d, 
and the fixed point value 
16 j „ f d\ d-2 



(107) 



fn*e. (108) 



In one and two space dimensions, we obtain A* = 32/3 
and A* = 32\/2/3, respectively; in three space dimen- 
sions, we recover the free theory on large scales. The 
nontrivial dependence of A* on the space dimensionality 
d, which is deeply related to the dissipative smoothing 
mechanism, is shown in Fig. |3J Note that the occurrence 
of q A in the denominator of the integrand in Eq. (|106l) is a 
direct consequence of using the Laplacian as the simplest 
and most natural scalar differential operator in defining 
the coupling operators in Eq. f|5 1 [) . It is quite remarkable 
that the detailed form of this integral leads to a vanish- 
ing A* in three space dimensions and to the well-known 
leading-order coefficient in the e expansion of A*. Equa- 
tions (|107p and (|108|) moreover lead to an explicit result 
for the /3 function (f3"0"l) which agrees with the famous re- 
sult for the ip 4 theory as, for example, given in Eq. (11.17) 
of [H or in Eq. (18.5.7) of [39|. 



V. REMARKS ON GENERALIZATIONS 

Our development of the quantum theory of a scalar 
field in d space dimensions has been based on pairs of 
adjoint operators a\, and a*, creating and annihilating 
field quanta of momentum k £ K d . We here consider 
generalizations based on on collections of adjoint opera- 
tors and creating and annihilating field quanta of 



momentum fc, where J is an additional index labeling the 
different types of quanta and the Fock space needs to be 
extended accordingly. For vector instead of scalar fields, 
as the simplest example, J labels the different spatial 
components. 

For bosonic fields, all creation operators commute 
among each other, and so do the annihilation operators. 
The only nontrivial commutation relations are 



} = 5j,pS(k-k'). 



(109) 



In terms of creation and annihilation operators, the 
Hamiltonian of the free field theory in the generalized 
setting is assumed to be of the form 



ff (0) = / uj k ai f aid d k 



(110) 



where a summation over the index J is implied by the 
Einstein summation convention and cok is a non-negative 
real-valued function with the symmetry property uj-k = 
LUk ■ The standard form is given by the relativistic energy- 
momentum relationship (|44p . with a mass parameter m. 
For massless fields, we simply have uik = \k\. The fact 
that we assume ojk to be independent of J (or that it 
depends in a restricted manner on J) expresses an un- 
derlying symmetry of the theory. In the same spirit, 
we can implement symmetry properties in the friction 
mechanism (however, Lorentz invariance would still be 
violated). 



A. Pure Yang-Mills fields 

For Yang-Mills fields, the label J for the field quanta 
consists of a spatial index and an additional label asso- 
ciated with the infinitesimal generators of an underlying 
Lie group, where we use the letters a, b, c for this kind of 
index [assuming 3 values for SU(2) corresponding to the 
W + , W~, and Z° bosons mediating weak interactions, 8 
values for SU(3) corresponding to the gluon "color octet" 
mediating strong interactions, and N 2 — 1 values for gen- 
eral SU(N)]. More precisely, for Yang-Mills fields, the 
components are associated with the base vectors T° of 
the Lie algebra of the underlying linear Lie group (we 
here consider matrix groups only [40l]). We assume the 
orthonormality conditions 



2tr(T a T b ) = S ab 



(111) 



The most famous example are the three Pauli matrices as 
base vectors T a of the Lie algebra of SU(2), which need 
to be divided by 2 to satisfy Eq. (|111[) . The base vectors 
T a are traceless matrices. 

In Yang-Mills theories, we deal with vector fields. The 
index J hence is of the form J = (j a) where j is a spa- 
tial index and a labels the base vectors of the Lie algebra. 
We use the amputated symbol j instead of j to emphasize 
that this spatial index takes only d — 1 instead of d val- 
ues. It is typical for gauge theories that one component 
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of a vector can be eliminated by gauge transformations 
so that vectors are reduced to d — 1 polarizations, here 
labeled by the amputated symbol j. For the field quanta 
this implies that the vector bosons mediating weak and 
strong interactions in d = 3 space dimensions come in 
two polarization states. 

The spatial components of the physical fields and their 
conjugate momenta are obtained as generalizations of the 
Fourier transforms (l40t and (14T1). 



— ik-: 



'^vWiv^'r - k ) 



and 



(112) 



n (ja) (a 



(113) 



where summations over j, Z are implied and the two po- 



larization vectors ej^ and typically depend on fc, 
yet in a symmetric manner. The polarization vectors 
express certain gauge conditions. We here choose the 
particularly convenient axial gauge of a vanishing last 
field component, A( da } — 0, because it avoids some sub- 
tle complications known for other gauges (see remarks 
on p. 15 of |39j and the added note in the conclusions 
of [4lj). For the axial gauge in d = 3 dimensions, we 
concretely choose [12] 



(2) 



J 1 ) 



k 3 v/fcf+fcf 



(114) 



and 



( 2 ) _ sgn(fc 3 ) 



-k 2 
A'i 




(115) 



These vectors ' and e^. ' are constructed such that 



and e^ 2) 

one obtains canonical commutation relations for the fields 

A (la) A (2a) and n ( lo ) n (2a) j 



[AM(x),nW(^] = i^j J i{(j 



(116) 



To obtain this generalization of Eq. (|42|) for arbitrary d, 
we postulate the property 



for j,Ze (117) 



j, but a further label a for the infinitesimal generator 
of the group occurs, and a so far unmotivated projector 
in momentum space has been introduced in Eq. (|113|) . 
Much deeper use of the properties of gauge theories must, 
of course, be made in the formulation of the Hamiltonian 
for the interactions between the various kinds of field 
quanta. As a first step, we need to look at the following 
constraint arising for Yang-Mills fields, 



djll^ + g f abc A^Xl { ^ = 0, 



(118) 



where g is the coupling constant and f abc are the struc- 
ture constants of the underlying Lie group, say SU(n). 
For the free theory with g = 0, the projector in Eq. (|1 13[) 
implies that the constraint (II 18[) is fulfilled so that its 
occurrence is now motivated. In the presence of inter- 
actions, we modify only the last component lT da ) such 
that Eq. (fTT5|> is fulfilled. The components 11^ for 
j < d — 1, which appear in the canonical commutation 
relations (|116[) . remain unchanged. 

The constraint (II 1 8[) . which fixes Il( da \ is a result 
of the field equations for the Yang-Mills theory. It de- 
pends on the strength of the interaction inherited from 
the Hamiltonian. This is very different from a gauge con- 
dition, say from fixing A^ da ^ = 0. The gauge condition 
is an example of a primary constraint, whereas Eq. (|118p 
is an example of a secondary constraint resulting from 
the consistency of the primary constraints with the evo- 
lution equations. Even more important is the distinc- 
tion between first and second class constraints. Second 
class constraints can be treated in a consistent canoni- 
cal quantization procedure based on the Dirac brackets 
of classical mechanics or field theory (see, for example, 
[HI or Sec. 7.6 of [3(3]). First class constraints are typ- 
ically related to unphysical degrees of freedom and can 
be taken into account by arbitrarily choosing a particu- 
lar gauge. Unfortunately, there is no guarantee for the 
gauge invariance of the final results obtained by canoni- 
cal quantization in different gauges. A loss of symmetry 
seems to be the unavoidable price to pay for a simple 
canonical quantization 44] . 

Within the axial gauge, the full Hamiltonian can be 
expressed nicely in terms of the spatial components of 
the fields A^ a ^ and IlW (see, for example, Eq. (15.4.10) 

of (ai) 



h = J in (3 ' a ) (x)nW°J (a) + ( X )F^ (x) 



d 3 x, 
(119) 



with the tensor fields 

p(ij<>) — d i A^ — d j + gf abc A^A^ c \ (120) 



for the polarization vectors of the axial gauge with e^ j = 
0, where summation over j but no integration over fc is 
implied in Eq. ([TT7j) . 

So far, we have made little use of the fact that we 
are dealing with a gauge theory: After fixing the gauge, 
J = (j a) has only d — 1 spatial components labeled by 



With the Fourier transforms (|112p , (|113[) , the gauge con- 
dition A( da ) = 0, and the modification of n^ da ' required 
to fulfill the constraint (|118[) . the Hamiltonian (|119l) can 
be rewritten in terms of the creation and annihilation 
operators a'^} and a'^. Except for an irrelevant constant 
contribution, the free Hamiltonian for g = is of 
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the form (jllOp provided that u>k — \k\, that is, Yang- 
Mills fields must be massless. The interaction Hamil- 
tonian contains first- and second-order terms in g and a 
few field operators only, which is an extremely convenient 
feature of the axial gauge, in particular, for perturbation 
theory. Note that the form of the Hamiltonian describing 
the interactions depends on the gauge. It even happens 
that the interactions are described by a simple polyno- 
mial Hamiltonian in one gauge and by a nonpolynomial 
Hamiltonian in another gauge [4l|, |42l EH . 

With the described steps to introduce fields and Hamil- 
tonian, it is possible to use the thermodynamic coarse- 
graining approach to the quantization of pure Yang-Mills 
fields. The introduction of friction simply happens by 
adding a label J (to be summed over) to each of the cre- 
ation and annihilation operators in Eq. (|54l) . Also the 
coupling to matter would be straightforward. On the 
one hand, the loss of gauge invariance in the proposed 
procedure admittedly is a high price to pay. The results 
for gauge-invariant correlation functions, such as Wilson 
loops, need to be compared to conventional calculations 
and the ft function needs to be evaluated. On the other 
hand, a major benefit may be the simultaneous solution 
of all ultraviolet and infrared problems, where the latter 
are known to be particularly subtle for the present ap- 
proach to Yang-Mills theories because axial gauges lead 
to additional singularities in the propagator so that naive 
principal- value recipes do not work [46145 Of and inconsis- 
tencies in the renormalization-group flow in perturbation 
theory might arise EH . 

B. Gravitation 

In the ip 4 and Yang-Mills theories, the friction mech- 
anism provides a dynamic cutoff on the length scale £ 
introduced in Eq. (|87|) . This smoothing length scale has 
no physical significance. We are dealing with an effective 
field theory valid only at length scales large compared to 
£. Physical predictions are obtained after removing the 
cutoff from the results of detailed calculations, that is, 
in the limit I — > 0. In the description of gravity, a quite 
different scenario should be expected. The smoothing 
length scale £ should be of the order of the Planck length 
and the dissipative friction mechanism should express a 
physical smearing of space and time at short scales, a 
fundamental effect that precludes further resolution. We 
hence expect an expression for the friction parameter 7 
of the form 



where G is the gravitational constant and a numerical 
prefactor still needs to be determined (we here show the 
occurrences of the reduced Planck constant and the speed 
of light explicitly). The determination of the prefactor 
could be based on a discussion of the entropy produc- 
tion rate as associated with black holes 153 55] or more 



general gravitational fields (for a thermodynamical dis- 
cussion of the entropy of classical gravitational fields see 
[56|). However, such wild speculations should be post- 
poned until we will have succeeded in formulating a con- 
vincing covariant friction mechanism. 

The above discussion implies that the role of 
renormalization-group theory is changed in quantum 
gravity. For the ip 4 and Yang-Mills theories, the ft func- 
tion characterizes the running coupling constant of a min- 
imal model that is appropriate to predict only the large- 
scale properties of its universality class. If the small reg- 
ularizing length scale used for the minimal model goes to 
zero, we obtain the usual rules for handling the singular- 
ities of a field theory. In gravity, on the other hand, the 
small length scale £ has a physical meaning. We need a 
proper physical model on the scale t, not just a minimal 
model. No limit is required, no singularities occur. As 
gravity can be observed over an enormous range of length 
scales, the coupling constant at the Planck scale must be 
at its critical value. The Planck scale, the gravitational 
constant, and the friction parameter hence contain equiv- 
alent information, as is clear from the definition of the 
Planck scale and Eq. (|121[) . For gravity, renormalization 
is not related to fundamental aspects of constructing a 
field theory in terms of minimal models but belongs only 
to the realm of perturbation theory. When we pass form 
the Planck scale to larger scales, the original physical 
model is degraded to a minimal model for which, on the 
scales of the order of the scales of physical interest, we 
can apply perturbation theory with the hope to get useful 
results. 

For the detailed mathematical formulation of the dis- 
sipative quantum theory of gravitation one might try to 
follow the procedure sketched for Yang-Mills theories. 
However, we are then faced with the same problems as in 
the canonical quantization of gravity. As a consequence 
of their complicated nonpolynomial functional form, the 
handling of the constraints resulting from gauge invari- 
ance is much more difficult for gravity. In particular, 
problems arise for the construction of the momentum 
representation. Ashtekar's famous idea of enlarging the 
gravitational phase space to incorporate spinors and to 
represent the spatial components of a metric in the de- 
composed form gij = —<JiA B PjB A (with summations over 
the SU(2) spinor indices A and B) in order to obtain sim- 
pler constraints [57], l58l (see, for example, also Sec. 3.2 
of [59| and Sec. 4.3 of |6fJ) may well be the key to obtain 
the dissipative quantum theory of gravitation. 



VI. SUMMARY AND CONCLUSIONS 

Within the framework provided by the thermodynamic 
quantum master equation, we propose to include irre- 
versible terms into quantum field theory and to consider 
the limit of weak friction and low temperatures. These 
irreversible terms account for the fact that we need to 
eliminate degrees of freedom below certain length and 
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time scales in order to be able to make use of the field 
idealization. The smoothing length scale resulting from 
the dissipative friction mechanism provides an alterna- 
tive regularization scheme where an external cutoff is re- 
placed by a dynamic process. However, the change from 
Hamiltonian to dissipative equations is much more than 
just another regularization scheme because it comes with 
a number of important implications. The vacuum state 
is no longer given by the ground state, but rather by a 
canonical density matrix characterized by a temperature 
T e . The Hamiltonian has a double role because it oc- 
curs both in the canonical density matrix and serves as 
a generator for the reversible contribution to dynamics. 
Entropy occurs naturally; it may be irrelevant in effec- 
tive field theories because the entropy production rate is 
negligibly small, but it must play an important role in a 
full quantum theory of gravity. 

In applying the proposed ideas to ip 4 theory in d space 
dimensions, we have have elaborated explicitly that a 
perturbation theory can be constructed with guidance 
from a detailed-balance principle, without any need to go 
through Dyson's U matrix to obtain meaningful results 
(see the discussion in Sec. 17.1 of 28]). All regularization 
is consistently provided by the irreversible contribution 
to time evolution that rapidly damps the local degrees of 
freedom. Although the irreversible friction mechanism is 
implemented in a non-covariant manner, the final physi- 
cal predictions of second-order perturbation theory in the 
limit of vanishing friction coincide with the well-known 
manifestly covariant results. A covariant formulation of 
the friction mechanism would be desirable but, in view 
of the second-order derivatives occurring in Eq. (|5ip , it 
presumably requires the introduction of additional fields 
corresponding to spatial derivatives 12 1. 

In the present approach, renormalization-group theory 
is used as a tool to refine perturbation theory. It requires 
a certain characteristic occurrence of the small smoothing 
length scale in perturbation expansions, which we have 
verified explicitly to second order. As a result, we repro- 
duce the well-known expression for the /? function for the 
running coupling constant. The fixed-point value of the 
dimensionless coupling constant is found to depend in a 
nontrivial way on space dimensionality. 

Although we focus on perturbation expansions, the ap- 
proach to quantum field theory proposed in the present 
paper is by no means restricted to perturbation theory. 
The fundamental quantum master equation (|54|) could 
also be treated by non-perturbative methods, including 
stochastic simulation techniques. Not even the lineariza- 
tion of the master equation around equilibrium would 
be necessary if one liked to benefit from the full struc- 
ture of the thermodynamic framework, for example, for 
fundamental or constructive developments. The underly- 
ing equation for the evolution of dissipative systems has 
a deep geometric structure and hence provides a sound 
basis for the description of dynamic quantum systems, 
including a Lyapunov function (entropy). The formu- 
lation is done according to the principles of nonequilib- 



rium thermodynamics as the proper framework for sys- 
tems with eliminated degrees of freedom. 

For proving the new approach to be useful, one should 
establish its applicability to a wide range of quantum 
field theories. We have discussed how dissipative smooth- 
ing can be applied to pure Yang-Mills fields where we 
propose to do all calculations in the axial gauge. As a 
crucial test, the proposed dissipative friction mechanism 
should resolve the subtle infrared problems occurring for 
Yang-Mills fields in the axial gauge by providing the clear 
causality properties of irreversible equations. For the cal- 
culation of Wilson loops, the proper definition of multi- 
time correlations needs to be understood, including the 
proper implementation of time-reversal symmetry, de- 
tailed balance, and their preservation under perturbation 
theory. For quantum gravity, the frictional smoothing 
mechanism would need to be elevated from a regulariza- 
tion scheme to a physical smearing of space and time. 

As a first step, we have restricted ourselves to second- 
order perturbation expansions in ip 4 theory. Of course, 
the dissipative approach to quantum field theory should 
also be validated for higher orders of perturbation the- 
ory and for more complicated field theories. One might 
want to attempt to establish a structure of perturbation 
expansions consistent with renormalization-group theory 
by analyzing Feynman diagrams. However, in view of 
the robust geometric structure of the thermodynamic ap- 
proach, a deeper discussion of the limits of vanishing fric- 
tion and temperature seems to be the more relevant issue. 
To simplify the practical calculation of perturbation ex- 
pansions, one needs to formulate rules to identify those 
Feynman diagrams that do not contribute in the limits 
of vanishing friction and temperature. Many details of 
dissipative quantum field theory still need to be clarified, 
many possibilities still need to be explored. 
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Appendix A: Some useful integrals 

We here collect a number of useful integrals that occur 
in the perturbation expansion of ip 4 theory. To facili- 
tate a direct comparison, we introduce the integrals 7i, 
I2, I3 in exactly the same way as in Sec. 11.5.1 of [34j . 
These integrals depend only on the dimensionality d. As 
before, we use d for the dimension of space, D = d + 1 
for the dimension of spacetime, and e = 3 — d = A — D. 
The D = d + 1 components of a spacetime vector k are 
given by k — (k, fc), where k is real and k has d com- 
ponents. Connecting integrals in d and D dimensions, 
where we consider both real and imaginary time, is the 
key to revealing Lorentz invariance. 
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The simplest integral is defined by (see Eq. (11.39) of 

34]) 
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Note that the integral I\ may be ill-defined in certain 
integer dimensions. If k corresponds to the real fre- 
quency, one needs to be careful with the proper handling 
of the poles which are located on the integration line; 
then, the is with infinitesimal positive e shifts the poles 
to avoid ambiguities. If K corresponds to the imaginary 
part of the frequency, the integration path is well sep- 
arated from the poles and the is prescription becomes 
irrelevant. After properly closing the contour for the k 
integration, Cauchy's integral theorem of complex anal- 
ysis allows us to pass from the D-dimensional integrals 
to the d-dimensional representation. 

By means of the general formula for reducing d- 
dimensional spherically symmetric integrals to one- 
dimensional integrals, 



f(k)d d k = 
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we obtain the explicit result 
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The following two useful integrals are closely related to 
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As a next step, we consider the dimensionless integral 

™6-2D 
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where, again, certain integer values of D should be ex- 
cluded. This integral has been discussed in great de- 
tail on pp. 282-284 of 34]. According to Eqs. (11.38), 



(11.40), and (11.41) of (34|, the useful integrals I 2 and I 3 
can be introduced by expanding Ik(u> 2 ) in terms of ui 2 , 



= h-h 



k=0 



0( w 4 ). 
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By deforming the integration path from imaginary to real 
frequencies and carefully stating by the is prescription 
how poles on the new integration line have to be cir- 
cumvented, we arrive at the Lorentz invariant expression 
occurring, for example, in truncated form for D — 4 on 
p. 174 of [H, 
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Again, we can pass from D-dimensional to ci-dimensional 
integrals. After replacing the factor 5(k\ + n 2 + k 3 + 
k) contained in Eq. (|A8|) by its one-dimensional Fourier 
transform, the integrations over K\, k 2 , and n 3 can be 
performed. After a further integration of an exponential, 
we obtain 



,4-2d 



d d hd d k 2 d d k 3 
(2n) 2d J 4uj kl uj k2 uJk 3 

Wfcj + Wfe 2 + LJk 3 



5(ki + k 2 + fe 3 + fe) 
(A9) 



(uj kl + tu k2 + uj k3 - is) 2 - LU 2 ' 

For the actual evaluation of Ik(u! 2 ), the form (|A6I) of- 
fers the most convenient starting point. Note that each 
of the factors in the denominator of Eq. (|A6j) contains 
a positive-definite quadratic form Q of the integration 
variables. After using the identity 1/Q = / °° e~ Qz dz 
for each of the three factors in the denominator, the 
two Gaussian D-dimensional integrations remaining af- 
ter making use of the 8 function can be carried out easily 
and only three one-dimensional integrations remain to be 
done (see Sec. 9.6 of [34]]). 

In comparing our results to those of [H|, for example, 
the formula for h given in Eq. (|A3I) to (11.44) of [34$ . 
the following identities are useful (see 8.334 of j38|): 



Y(l-x)Y{x) 



sin nx 



(A10) 



(All) 



Appendix B: Some useful examples of Ta 

Among the normal ordered products of up to three 
creation/annihilation operators, Ta vanishes for all op- 
erators except for a^. a,k 2 , a]^ afe 2 afe 3 , and a£ a^afe 3 - For 
the product of two operators, we have 

r <ah2 = wM ^6( kl -k 2 ). (Bi) 
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For the product of three creation/annihilation operators, 
we find 



r t 

a kl a k 2 a h 3 



and 



|^(fci-fe 2 ) W(uj kl ,-uj kl - uj ka ) 
+ W(-uj kl ,uj kl -uj k3 ) a k3 
+ S(k 1 -k 3 ) W(uj kl ,-uj kl - uj k2 ) 
+ W(-uj kl ,uj kl -Wfe 2 ) a fe2 |, (B2) 



^S{k 2 -k 3 ) W(u ks ,u kl - w fe3 ) 

+ W(-uj k3 ,u kl + uj k3 ) a' ki 

+ S(k 1 -k 3 ) W(uj k3 ,uj k2 ~ uj k3 ) 

+ W{-uj k3 , u k2 + uj k3 )] a{ 2 } . (B3) 



3 - 1 



These results are obtained by straightforward evaluation 
ofEq. GSJ). 



J 



Appendix C: Some useful commutators 

To facilitate the calculation of averages occurring in perturbation theory, we calculate some basic commutators 
involving the Hamiltonian. They are obtained in a straightforward manner from the fundamental commutation 
relations (j42|) and the from (|45|) of the ip 4 Hamiltonian. The following single commutators are used in evaluating the 
averages in Eqs. (J78]) and (|T03| : 



[a h ,HW] = ¥-(a k + al k 

U) k \ 



A 1 



1 f d d kid d k 2 d d k 3 



<5(fci + k 2 + k 3 + k) 



x a 



24 (27r) d y/u^ J ^/Q kl Cj k2 ui k3 
kl a-k 2 a- k3 + Sa^a-^a-^ +ial i a { k2 a- k3 + 4i a L a L)< 



(CI) 



and 



d d k 1 d d k 2 d d k 3 



<5(fci + k 2 + k 3 - k) 



24 {2it) d J y/Q kl uj k2 Lj k 
x (a- kl a- k2 a- ka + ia^a^a^ + 3a* ± al 2 a-k 3 + <4.! a L a fe. 



(C2) 



The following type of identities for single commutators allows us the separation of the two contributions represented 
by the Feynman diagrams in Fig.[^a) and (b): 



and 



[a kl a k2 ,H {1 '>} = [a fcl ,iJ (1) ]a fc2 + [a k2 , H (1) ] a kl + [a kl , [a k2 , H (1) }}, 



(C3) 



(C4) 



We also provide a number of double commutators. The following double commutators are used in Eqs. (|77p and 
(H02]): 



[a kl ,[a k2 ,H^] 



A 1 



1 



f d d k[d d k 2 :„,.,, , s ( , t 

/ 5 <n«i +k 2 + kx + k 2 ) a_ fc ' a_ fc < + 2a' , a_ fcf 



8 (2^) d ^ k ~^ k ~ 2 

^As(ki + k 2 ), 
u kl 



y/^k'^k 



+ at, at, 

n, 1 ft- 2 



(C5) 



[[H 



(i) „t l „t i = 



A 1 



(27r) rf y/Q kl Cb k2 J ^wt.;^ 
Az 

— o(fci + fc 2 ), 



5(fci + fc 2 - ki - k 2 ) (a-k^a-k^ + 2a* h ,ji- k ' 2 + a£,a£. 



(C6) 
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and 



[a kl ,[H^,a{ ]] = ± 1 [ d ^p_ 6 (k[ +k' 2 + k 1 - fc 2 ) (a_ K a_ K + 2o^ .a_ K + a{, a{^\ 



(C7) 



The following type of identities for double commutators facilitates the evaluation of the average in Eq. (|100p : 

[ak 1 a.k 2 ,[ H{ ' 1 \ a \' 1 a \'J\ = 

S(k x - k[) a k2 [H^,al 2 ] + 6(ki - k' 2 ) a k2 a^J + S(k 2 - k[) a kl [H (1 \ a£, ] + S(k 2 - k' 2 ) a kl a£,] 

+4; [ofc 3) [-ff (1) ,a^]]afei + 4 2 t afc ^ 4'J] a fei + 4; Ki > [# (1) > 4;J] a fc2 + 4 2 [ahi > [# (1) > 4'J] a &2 

+4' [afcu [dfc 2 , [if , a+ , ]]] + o^, [a kl , [a fc2 , [if « , 4- ]]] + [a fcl , [[ff (1) , a{, } , a{, }} a k2 + [a k2 , [[fT« , a+, ], 4 , ]] a fel 

+ [a fel ,[a fe2 ,[[ff«,at,],a t fe ,]]] ! (C8) 

and 

[[a kl a k2 ,H { %al ti al, 2 ] = 

S(kx - k[) [a k2 , H «]a£, + 5(fer - feQ [a fe2 , Wfcj^ + 5(fc 2 - fc£) [a fcl , «]^, + <5(fc 2 - fe 2 ) [a kl , < x )]o+, 
+4; [afc 2 , [i? (1) ,a^]] a fcl + a£, [a fc2 , a^J] a fcl + 4; [a fel , a^]] a k2 + 4 2 [a fel , a£, ]] a fc2 

+0^/ [ofcn[afe a ,[J (1) ,ay]] + 4 2 t afc i' K 2 > > a fe'J]] + [ a fei 3 [[- ff(1) : a fc' 1 ]'«l 2 ]] a fc 2 + i a k 2 , [[-H" (1) ,ai'J' a l 2 ]] afe i 
+[o fcl ,[o fca ,[[FW,oi,],oL]]]. (C9) 



Appendix D: Exact second-order perturbation results for propagator 

By evaluating the average ((77)) by means of the first order expansion (|2T|) of averages, the double commutator (|C7I) . 
and Wick's theorem (l57l). we obtain the exact result 



Xi. — A 



J + 2z 
2Cb k 



A2 L_ /" A 

^2 



16w fc (2^) d 7 c; 



w(2w g ) 



lO(-2Wg) 



with 



J = 



(l_ e -/3w a )2 ( e/ 3c, e -^«) (e^-1) 2 



(27r) d 7 2Q q (eP u < - 1) ' 



(Dl) 



(D2) 



Note that, in the low-temperature limit f3 — > oo, the integral J becomes exponentially small. In a similar way, 
the average ([75)1 can be evaluated by means of the commutators (jCll) and (|C2[) . the formula (|75|) with the explicit 
expressions of Appendix IB1 and Wick's theorem (f57|) . The complete result is 



Y k (u) = -X 



2 z(J + 2z) 



J + 2z 1 

A 



1 



1 



d d q ( 1 



u\ (2Tr) d J Q q (e( lu *-l)\w + w k -i'y l A u - u k - n a t t „ 



, 2 J + 2z 1 
+ A 



8^ (2^) d 



?7 g [W^tJg, + a; fc ) + W^a^, -ui q - ui k )] d d q 
Cj q {e^" - 1) 



1 



1 



w + uj k - z 7a t w + w fe - z 7fc 



d d kid d k 2 d d k 3 



5(fci + fc 2 + fe 3 + fe) 



96^(2^)^7 Cj kl Cj k2 uj H (l_ e -' 3 ^i)(l-e-' 3 ^2)(l 



-/9wfc 3 
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\ _ g — /S(wfc 1 +Wfc 2 +Wfc 3 ) 



2 _ e -P(uik 1 +^k 2 +^k 3 ) 



UJ 



+ 3 



+ UJ kl + UJ k2 + U>k 3 - ha-^a. 
P -Pu k . _ -f3(u> k +u> k ) 



UJ - UJ kl 



Uk 2 



uj — uj k , + L0u o + wto — ij t 



-3 



-P^k 3 _ e -/3(w fcl +w fc2 ) 



w - w fcl - w fc2 + Wfc3 - n a t t a _ 

fcl >B2 



(D3) 
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